N94- 25130 


Thermophysical Properties of Simple Liquid Metals: A Brief 

Review of Theory 

David Stroud 
Department of Physics 
The Ohio State University 

In this paper, we review the current theory of the thermophysical properties of simple liq- 
uid metals. The emphasis is on thermodynamic properties, but we also briefly discuss 
the nonequilibrium properties of liquid metals. We begin by defining a “simple liquid 
metal” as one in which the valence electrons interact only weakly with the ionic cores, so 
that the interaction can be treated by perturbation theory. We then write down the equilib- 
rium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the 
electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of 
electrons and ions. Since the electron-ion interaction can be treated by perturbation, the 
electronic part contributes in two ways to the Helmholtz free energy: it gives a density- 
dependent term which is independent of the arrangement of ions, and it acts to screen 
the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density- 
dependent, in general. After sketching the form of a typical pair potential, we briefly enu- 
merate some methods for calculating the ionic distribution function and hence the Helm- 
holtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, 
and thermodynamic perturbation theory. 

The final result is a general expression for the Helmholtz free energy of the liquid metal. 
It can be used to calculate a wide range of thermodynamic properties of simple metal liq- 
uids, which we enumerate. They include not only a range of thermodynamic coefficients 
of both metals and alloys, but also many aspects of the phase diagram, including freez- 
ing curves of pure elements and phase diagrams of liquid alloys (including liquidus and 
solidus curves). We briefly mention some key discoveries resulting from previous appli- 
cations of this method, and point out that the same methods work for other materials not 
normally considered to be liquid metals (such as colloidal suspensions, in which the sus- 
pended microspheres behave like ions screened by the salt solution in which they are 
suspended. 

We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties 
which can be treated by an extension of these methods. These include electrical resistiv- 
ity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffu- 
sion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly 
mention two methods by which the theory might be extended to non-simple liquid metals: 
these are empirical techniques (i.e., empirical two- and three-body potentials), and 
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numerical many-body approaches. Both may be potentially applicable to extremely com- 
plex systems, such as nonstoichiometric liquid semiconductor alloys. 


What is a “Simple Liquid Metal”? 

A simple liquid metal is one in which electrons can be 
unambiguously divided into two categories - “valence” 
electrons and “core” electrons - such that the interac- 
tion between the valence electrons and the ionic core is 
sufficiently weak to be treated by low-order perturbation 
theory. Examples: alkali metals, alkaline earths, several 
polyvalent metals (such as Al, In). The transition and 
noble metals are not simple metals. 
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Hamiltonian of a simple liquid metal is 

1 Z 2 e 2 1 e 2 

H ~ 2 R^R. |R - R'l + 2 |r, - Fjl 


+ £ V'p.dr, - R|) 

r i,R 


Pi -Pr 

i2m r2M 


• R = ionic position (ions of valence Z) 

• P R = ionic momentum 

• T{ = electronic position 

• pi = electronic momentum 

• V p3 = electron-ion interaction ( “pseudopotential” ) 

F — Helmholtz free energy = kj?T In Z 

Z = ^N / n R dRdP R Tr e exp( - ~ ) 

Here h = Planck’s constant, “Tr e ” denotes a trace over 
electronic coordinates. The ions axe taken to be clas- 
sical particles, while the valence electrons are assumed 
quantum-mechanical. 
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Rewrite partition function Z as: 

Z = J^J n R dRdP R e- l3H ‘ i [Tr e e-K H ‘' +H ^} 

• H a = first and fifth terms in Hamiltonian on previous 
page 

• H ee = second and fourth terms 

• H e i = third term 

• P = l/k B T. 

We can now rewrite the electronic trace as 
Tr e [...\ = exp(-/?F'({R}), 

where F' e is the free energy of an electron gas in the “ex- 
ternal potential” produced by N ions located at (Ri,..., 
R n)- 
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In turn, 

F' e = F eg + Ebs 

• F eg = free energy of a uniform interacting electron 
gas in a uniform positive background. 

• Ebs = “band-structure energy” = extra piece in free 
energy due to electron-ion interaction. In spite of its 
name, the band-structure energy does not imply a 
periodic lattice. 

If V p3 is weak, we can do perturbation theory to estimate 
E55. Band-structure energy in weak-pseudopotential limit 
is 

Ebs » \ T, Xe(q)|V r p «(q)p(q)| 2 
2 q ^0 

*•<■» - 5? t® ■ - ■) 

p(q) = Ee- iqR , 

R 

where e(q) = dielectric function of interacting electron 
gas. (q = 0 term is not small and treated separately.) 
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We can now express the partition function in terms of 
ionic coordinates alone: 

Z = ^n! HRdRrfPR exp{-/ 3 [W e// ({R}, n e ) + e||]}, 

where W e /f is an effective potential energy, which has 
the form 

1 Z 2 e 2 

W eff ({K,n e }) = - E m — . + E bs ({R} + F eg (n e ). 

* R^R' |xt — K. I 

Note that E55 is a density-dependent, structure-independent 
part, plus a sum of two-body contributions (if V ps is 
weak). This leads to effective pairwise interionic interactions . 
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Schematic of a Typical Pair Potential 



nearest-neighbor separation 

• R h s w effective (temperature-dependent) hard-sphere 
diameter ~ distance of closest approach at tempera- 
ture T. 

• Rjon — ionic diameter (R lon < R hs)- 
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Given the pair potential, the ionic distribution function 
can be calculated by several methods. For example: 

• Monte Carlo simulation: a statistical sampling of 
configuration space, with weight proportional to the 
Boltzmann factor; 

• Molecular dynamics: solve Newton’s equation of mo- 
tion numerically at a given energy E (or temperature 

T); 

• Thermodynamic perturbation theory (approximate 
the liquid structure factor by that of a known liquid 
system (e. g., a liquid of hard spheres), then obtain 
“best” parameters (e. g., the hard-sphere diameter) 
variationally. This procedure gives a temperature 
and density-dependent effective hard-sphere diame- 
ter. 


Results 

F = Em + Ebs + E 0 — T 5 1<m 4- F eg + -NksT 

• E m = Madelung energy (bare ion-ion interaction en- 
ergy) 

• Ebs = band-structure energy 

• E 0 = “Hartree” energy (sum of long- wavelength terms 
in potential energy) 

• = ionic entropy 

• F eg = free energy of uniform, interacting electron gas. 
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Some Calculable Quantities 

• Free energy F = E - TS 

• Internal energy E 

• Entropy S 

• Specific Heat C y 

• Thermal expansion coefficient a = y (|^) 

• Isothermal compressibility kt = {jp) T 

• Liquid structure factor S(q) 

• Equation of state P(V, T) 

• Heat of formation (of elements) 

• Heat of mixing (of alloys) 

• Freezing curve (of pure elements) (i. e. line in P-T 
plane on which the liquid element freezes) 

• Properties on the freezing curve (A S, A V) 

• Phase diagrams of liquid alloys (phase separation 
curves in concentration-temperature plane, liquidus, 
solidus, eutectic points) 
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Some Comments on Specific Properties 


• Cy generally dominated by ions . (Valence electrons 
are in their ground state, or nearly so, at least for 
simple liquid metals.) 

• kt dominated by electrons . (Here, the density-dependent, 
structure-independent part of the free energy is cru- 
cial.) 

• Melting temperatures of pure elements reduced by 
screening (E bs )- (Ebs= 0 corresponds to the well- 
known “one-component plasma”.) The reduction is 
about 30% for the alkali metals, more than 80% for 
most polyvalents (Al, In, etc.). 

• Same theory for freezing also works for “colloidal 
crystals” (suspensions of charged polystyrene micro- 
spheres in a salt solution (Shih and Stroud, 1983). 
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Non- equilibrium (transport) properties of 
simple liquid metals 

• DC electrical resistivity. Reasonable agree- 
ment with weak-scattering (Ziman) formula 

p oc J 0 2kF k z \v;?(k)\ 2 s(k)dk 

where S(k) is the structure factor and V*™ (k) is the 
screened electron-ion pseudopotential. 

• Thermal conductivity. There are two contribu- 
tions: (i) electronic; (ii) ionic. Electronic contribu- 
tion dominates in simple liquid metals, and can be 
calculated from p using Wiedemann-Franz law (this 
is not true in poorly-conducting liquid semiconduc- 
tors). 

• Viscosity, Atomic Self-Diffusion Coefficient, 
Concentration Diffusion Coefficients in Al- 
loys. All these quantities are dominated by the 
screened pair potentials described above, and can be 
calculated from them by molecular dynamics simu- 
lations, or estimated from hard-sphere results (using 
variationally determined hard-sphere diameters). 
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Other Properties. 


• Surface Tension 7 . Reasonable results can be 
obtained from an “ionic density functional theory” 
which considers the free energy of an inhomogeneous 
liquid in terms of the free energy and response func- 
tions of a homogeneous liquid (Wood and Stroud, 
1983; Zeng and Stroud, 1986). Good results also for 
surface entropy — 

• Thermal Emissivity. This requires knowledge of 
reflectivity R(o>) (or equivalently, complex dielectric 
function e(w)) as a function of frequency in the in- 
frared. The Drude model is probably inadequate, 
even for simple metals; one should include band- 
structure effects at least in perturbation theory. 
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Extension to Non-Simple Metals (or to 
Liquid Semiconductors) 

There are several possibilities: 

• Use empirical (two-body and three-body) potential 
fitted to appropriate experiments. This is not very 
well justified from first principles, but may work in 
selected noble and transition metals, or liquid semi- 
conductors. 

• Use numerical many-body techniques. This tech- 
nique was pioneered by Car and Parrinello. It is very 
computer-intensive, but may be the best method for 
extremely complex systems, such as nonstoichiomet- 
ric liquid semiconductor alloys. With the advent of 
large computers, it may become feasible in the near 
future, both for thermodynamic and for transport 
properties. 

Conclusion 

Theory is rapidly improving! 


176 


